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ABSTRACT 

Observations of the high energy sky, particularly with the INTEGRAL satellite, have quadrupled 
the number of supergiant X-ray Binaries observed in the Galaxy, and revealed new populations of 
previously hidden High Mass X-ray Binaries (HMXBs), raising new questions about their formation 
and evolution. The number of detected HMXBs of different types is now high enough to allow us to 
carry out a statistical analysis of their distribution in the Milky Way. For the first time, we derive the 
distance and absorption of a sample of HMXBs using a Spectral Energy Distribution fitting procedure, 
and we examine the correlation with the distribution of Star Forming Complexes (SFCs) in the Galaxy. 
We show that HMXBs are clustered with SFCs with a typical cluster size of 0.3 ± 0.05 kpc and a 
characteristic distance between clusters of 1.7 ± 0.3 kpc. Furthermore, we present an investigation of 
the expected offset between the position of spiral arms and HMXBs, allowing us to constrain age and 
migration distance due to supernova kick for 13 sources. These new methods will allow us to assess 
the influence of the environment on these high energy objects with unprecedented reliability. 

Subject headings: Galaxy: stellar content - Stars: distances, early-type, evolution, formation - X-rays: 
binaries 



1. INTRODUCTION 

X-ray binaries are separated in two classes: Low Mass 
X-ray Binaries (LMXBs) and High Mass X-ray Binaries 
(HMXBs), differing by the mass of the companion star 
and the way accretion of matter occurs. HMXBs are 
binary systems composed of a compact object (neutron 
star or black hole candidate), accreting matter from a 
massive companion star, either a main sequence Be star 
or an evolved supergiant O or B star. Be stars are sur- 
rounded by a circumstellar decretion disc of gas, and 
accretion occurs periodically when the compact object 
passes through this disc of matter, whereas supergiant 
X-ray binaries are mostly wind-fed systems (see Chatyl 
Sonjfor a review about these different kinds of HMXBs) . 
Most of these source s are observed in the Galactic plane 
( Grimm et al.||200^ ) as expected since short-lived stellar 
systems do not have time to migrate far from their birth- 
place. 

Due to dedicated observations with RXTE and INTE- 
GRAL^ around 200 HMXBs are currently known in the 
Milky Way allowi ng us to focus on thei r distribution. Us- 
ing RXTE data, [Grimm et al] ( |2002D highlighted clear 
signatures of the spiral structure in the spatial d i stribu - 



tion of HMXBs. I n the same way, iDean et al. J2005), 



Lutovinov et al. 



(20051), Bodaghee et al.l (120071) and 



Ebdaghee et al. (2012 ) showed that HMXBs observed 
with INTEGRAL also seem to be associated with the 
spiral structure of the Galaxy. However, HMXB loca- 
tions, mostly derived from their X-ray luminosity, are 
not well constrained and highly uncertain due to direct 
accretion occurring in HMXBs, well below the Edding- 
ton limit. 



In order to overcome this caveat we present a comprehen- 
sive and novative approach allowing us to derive the loca- 
tions of a sample of HMXBs, and study their distribution 
in the Galaxy. We build the Spectral Energy Distribution 
(SED) of each HMXB and fit it with a blackbody model 
to compute the distance of each source. We study this 
distribution and the correlation with Star Forming Com- 
plexes (SFCs) observed in the Galaxy. Knowing the loca- 
tion of these sources, one can examine the composition of 
the environment at their birthplace. This st udy, whose 
preliminary results have been presented in Coleiro fc 
Chaty (2011), is necessary to better understand the for- 
mation and evolution of the whole population of HMXBs 
and primarily to state the role of the environment and 
binarity in the evolution of these high energy binary sys- 
tems as outlined in the last section of this paper. In 
Section 2 we explain how we derive HMXB distances, 
then, in Section 3 we show that HMXBs are correlated 
with Star Forming Complexes. In Section 4, we derive 
the expected offset between HMXBs and Galactic spiral 
arms before discussing some implications on HMXB for- 
mation and evolution. We also derive the age and kick 
migration distance for 13 sources. Finally, we conclude 
in Section m 



2. 



DERIVING THE HMXB LOCATION WITHIN OUR 
GALAXY 



To compute HMXB distances, we gathered a sample 
of HMXBs for which at least 4 optical and/or near 
infrared (NIR) magnitudes were known. For each source 
we build its SED and fit it with a blackbody model. 
This enabled us to evaluate the distance of the source 
along with its associated uncertainty. 
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There are currently more than 200 HMXB s detected 
in the Milky Way. Using the Liu et al. (2006) catalogue, 
updated with literatu r e, and the INTEGRAL source cat- 
alogue of Bird et al.' (120101), we retrieved the quiescent 
optical and NIR magnitudes (mostly from the 2MASS 
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point source catalog), the spectral type and the lumi- 
nosity class of each source from the literature. For each 
HMXB, four magnitudes are required for the fitting pro- 
cedure, since the condition N — n > 2 with N the num- 
ber of observed magnitudes and n the number of free 
parameters (in our study, n = 2, cf. section 2.1) needs 
to be met for statistics. Finally, 59 sources follow 
these requirements. 9 of them appear to be located in 
the Magellanic Clouds and four others are poorly known 
systems (IR counterparts are debated or physical param- 
eters were not found in the literature) so our final sample 
consists of 46 sources in our Galaxy which meet the pre- 
vious conditions (see Table [T]). 



2.1. SED fitting procedure 

Our fitting procedure is based on the Levenberg- 
Marquardt least-square algorithm implemented in 
Python (least sq routine from the scipy package). For 
each HMXB, we build the SED in optical and NIR from 
a maximum of 8 (U, B, V, R, I, J, H, Ks) and a minimum 
of 4 magnitudes points. This SED is then fit (see Figure 
IT] for an example) by a blackbody model, with the aim 
of using an uniform fitting method for the entire sample 
of HMXBs, given by the relation: 



A4 



X 10 



-0.4A;, 



(R/D)^ 



exp 



(1) 



where: A is the wavelength in /im; Fx the fiux density 
in Wm~^/im~^; h the Planck constant; c the speed of 
light; Aa the extinction at the wavelength A; R/D the 
stellar radius over distance ratio and ks the Boltzmann 
constant and T the temperature of the star. When the 
radius and the temperature of the companion star are 
available in the literature, we use these values (given in 
Table |2|. For the other sources, we derive the radius and 
the temperature of the companion star, which dominates 
the optical and NIR flux, from the spectral type and the 
luminosity class, by using tables of ,Vacca et al.^,1996 
Panagia||1973| [Martins fc PleSl[2QQ6ll5earle et al.||2QQ8 
thereafter noted PVMS. Since only four sources (includ- 
ing GX 301-2) with a highly ev olved mass star h ave been 
detected in the Milky Way (Mason et al. 2012), we can 
reasonably assume that for all the other systems studied 
here, the radius of the companion star is sensibly close 
to the one expected, depending on the luminosity class. 
When the spectral type is poorly known, we derive the 
mean distance value by considering the different possi- 
ble spectral types and we allocate adequate error on the 
distance of these systems. Finally, if no data is given 
" " ^r^.^ .^^ ^r^^^ "Searle et al 2008] 

is the case for lu- 
minosity classes II and IV) we computed intermediate 
R and T between luminosity classes I and III for class 
II and between class III and V for class IV respectively. 
This affects only 3 sources in our sample. For stars with 
peculiarities (N III and He II in emission mentioned by 
an 'f in the spectral type, broad absorption mentioned 
by an 'n', unspecified peculiarity mentioned by a 'p') we 
assume their radius and temperature to have the same 
values than "normal" stars. For Be stars, we take R and 



in Vacca 'eTari [19961 [Panagia[ [T973| [5 

about a certain spectral type (which i 



T values of normal B stars, and the circumstellar disc 
emission is taken into account through a possible uncer- 
tainty due to the expe cted infrared excess of this family 
of sources (see Section 2.2). Finally, we describe our es- 
timate of distance uncertainties due to errors on R and 
T of companion star in Section [2^ 



Moreover, two parameters are left free: the extinction 
in V band Ay and the ratio R/D, whereas the extinction 
A\ is derived at each wavelength from Cardelli et alj 



(1989) assuming Ry = 3.1 (corresponding to the mean 
value derived for a diffuse interstellar medium towards 
the Galactic plane 



Fitzpatrick[[1999 ). 

of potential change of the extinction 



The influence 

aw was explored. 

The results, presented in Section 2.2, do not show any 
substantial variation in distance determination due to 
extinction change. 

Knowing the radius R of the companion star, we then 
calculate the distance D in kpc. 



The 



X 



least square 



ction given by the formula 

r 1 2 

i,mode? (with Xi^ohs the observed 
flux value for the fllter z, Xi^^nodeh the predicted flux in 
the i^^ fllter derived from the black-body model and cr^, 
the flux error in the same fllter) is then minimized by 
the Levenberg-Marquardt algorithm. 

To check our results, we estimated the distance of each 
HMXB using the expression 



Dpc = IQ-^-^^"^^- 



-Av+5) 



(2) 

with mv the relative magnitude and My the absolute 
magnitude while Ay is given by the flt. This formula 
has the advantage of not depending on the companion 
star radius. However, it does depend on the a bsolute 



magnitude My, given in [Martins fc Plez[ pOOGb for O 
type s tars and in Morton ^ Adamsl ('1968* and[p; 



(1973) for B stars! 'I'he results obtained 



anagia 
with thii 

second method are consistent with those using the flrst 
approach. 

We give results in Table [l] and compare these results 
with previously published results in Appendix ??. We 
point out that median discrepancy in distance is ~ 17% 
and Ay is often very similar (median discrepancy of less 
than 7%). 



2.2. Uncertainties in the computed distance and 
extinction 

The magnitude uncertainties are retrieved from the 
literature. For sources for which no error is given, we 
use a systematic error of 0.1 magnitude, conservative of 
typical errors. The flux uncertainties are then derived 
from these magnitude errors. Nonetheless, we assume 
in the following that the spectral type given in the 
literature is the real spectral type of the companion star. 
We carried out simulations on a wide range of spectral 
types which enabled us to constrain the uncertainty that 
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NAME 


D (kpc) 


Av (mag) 


T (K) 


R (Rq) 


SpT 


lA 0535+262 


3.8 ± 0.33 


1.9 ± 0.26 


32930 


14.7 


09.7IIIe (09.5IIIe) 


lA 1118-615 


3.2 ± 1.4 


4.6 ± 1.9 


32930/34620 


14.7/8.50 


09.5IIIe/09.5Ve (09.5III/09.5V) 


IE 1145.1-6141 


10.5 ± 0.90 


5.2 ± 0.19 


18300 


51.0 


B2Iae (B2Ia) 


IH 1249-637 


0.63 ± 2.5 


1.8 ± 2.5 


30160 


14.8 


B0.5IIIe (B0.5III) 


IH 1555-552 


0.89 ± 0.093 


2.7 ± 0.6 


19500 


9.80 


B2IIIn (B2III) 


3A 0114+650 


6.5 ± 3.0 


4.0 ± 0.60 


24000 ± 3000 


37.0 ± 15.0 


Blla 


3A 0726-260 


5.0 ± 0.82 


2.5 ± 0.25 


38450/35900 


9.30/8.80 


08Ve/09Ve (08V/09V) 


3A 2206+543 


3.4 ± 0.35 


1.8 ± 0.60 


34620 


8.50 


09.5Vp (09.5V) 


4U 1700-377 


1.8 ± 0.15 


2.0 ± 0.15 


40210 


21.2 


06.5Iaf (06.5Ia) 


Cep X-4 


3.7 ± 0.52 


5.3 ± 1.4 


22600/20500 


6.17/5.62 


BlVe/B2Ve (B1V/B2V) 


Cyg X-1 


1.8 ± 0.56 


3.4 ± 0.18 


32000 


17.0 


09.7Iab 


EXO 0331+530 


6.9 ± 0.71 


6.0 ± 0.50 


37170 


9.00 


08.5Ve (08.5V) 


EXO 2030+375 


3.1 ± 0.38 


12 ± 1.4 


33340 


8.30 


BOVe (BOV) 


gam Cas 


0.17 ± 0.50 


1.2 ± 2.2 


25000-30000 


10.0 


B0.5IVe 


GRO J1008-57 


4.1 ± 0.59 


6.7 ± 1.1 


33340 


8.30 


BOe (BOV) 


GT 0236+610 


1.8 ± 0.20 


3.8 ± 0.65 


333340 


8.30 


BOVe (BOV) 


GX 301-2 


3.1 ± 0.64 


6.3 ± 0.14 


20400 


62.0 


Blla 


GX 304-1 


1.3 ± 0.10 


6.0 ± 0.17 


20500 


5.62 


B2Vne (B2V) 


Ginga 0834-430 


7.1 ± 4.2 


11 ± 2.2 


31540/33340/ 


14.7/8.30/ 


B0/2IIIe/Ve (B0/2III/V) 








19500/20500 


9.77/5.62 




H 0115+634 


5.3 ± 0.44 


6.4 ± 0.28 


33340 


8.30 


B0.2Ve (BOV) 


H 1145-619 


4.3 ± 0.52 


1.7 ± 0.55 


31540 


14.7 


B0.2IIIe (BOIII) 


H 1417-624 


7.0 ± 0.74 


6.1 ± 1.1 


22600 


6.17 


BlVe (BIV) 


H 1538-522 


6.2 ± 1.8 


6.4 ± 0.28 


28000 ± 2000 


17.0 ± 2.00 


BOIab 








31500 ± 1000 


17.0 ± 2.00 




IGR J00370+6122 


3.4 ± 1.2 


2.4 ± 0.19 


24300/30160 


21.75/14.8 


B0.5II/III 


IGR J01583+6713 


4.1 ± 0.63 


4.7 ± 0.32 


20000 


7.70 


B2IVe 


IGR J06074+2205 


4.5 ± 0.36 


3.3 ± 0.22 


32060 


8.00 


B0.5Ve (B0.5V) 


IGR J08408-4503 


3.4 ± 0.35 


1.7 ± 0.60 


34230 


23.8 


08.5Ib(f) (08.5I) 


IGR J11215-5952 


7.3 ± 0.68 


2.6 ± 0.60 


22000 


36.5 


Blla 


IGR J11305-6256 


3.6 ± 0.71 


1.0 ± 2.2 


31540 


14.7 


BOIIIe (BOIII) 


IGR J11435-6109 


9.8 ± 0.86 


5.7 ± 0.28 


32060 


8.00 


B0.5Ve (B0.5V) 


IGR J16465-4507 


12.7 ± 1.3 


5.0 ± 0.75 


23600 


33.1 


B0.5I 


IGR J17200-3116 


10.4 ± 3.6 


6.6 ± 4.5 


32060 


8.00 


B0.5Ve (B0.5V) 


IGR J18214-1318 


10.6 ± 5.0 


14 ± 3.1 


32740 


24.6 


09I 


IGR J18410-0535 


7.8 ± 0.74 


6.1 ± 0.65 


21700 


34.9 


Bllb 


IGR J18450-0435 


6.4 ± 0.76 


6.7 ± 0.49 


31240 


25.4 


09.5I 


KS 1947+300 


8.5 ± 2.3 


4.0 ± 0.49 


33340 


8.30 


BOVe (BOV) 


roti JDizo9-6o 


1.7 ± 0.56 


O.o ± U.7U 


Qonnn+2000 
32000_;l000 


9-00-1.5 


B2Ve 


RX J0440. 9+4431 


2.9 ± 0.37 


2.9 ± 0.25 


33340 


8.30 


B0.2Ve (BOV) 


RX J0812.4-3114 


8.6 ± 1.8 


2.3 ± 0.20 


28000 ± 2000 


10 ± 2.0 


B0.2IVe 


RX J1744.7-2713 


1.2 ± 0.46 


2.7 ± 0.57 


30160/32060 


14.8/8.00 


B0.5IIIe/Ve (B0.5III/V) 


SAX J1818.6-1703 


2.7 ± 0.28 


1.6 ± 0.80 


17000 


8.71 


B3III 


SAX J2103.5+4545 


8.0 ± 0.78 


4.2 ± 0.25 


33340 


8.30 


BOVe (BOV) 


Vela X-1 


2.2 ± 0.22 


2.2 ± 0.46 


24700 


33.8 


B0.5Iae (B0.5Ia) 


XTE J1855-026 


10.8 ± 1.0 


5.8 ± 0.90 


28100 


26.9 


BOIaep (BOIa) 


XTE J1946+274 


6.2 ± 3.0 


6.9 ± 0.74 


32430/22050/ 


11.5/8.10/ 


BO/l/IVe/Ve (BO/l/IV/V) 








33340/22600 


8.30/6.17 




X Per 


1.2 ± 0.16 


0.81 ± 0.22 


33340 


8.30 


BOVe (BOV) 



SpT Ref. 

Giangrande et al. 1 1980^ 
Janot-Pacheco et al. 1 1981. 
Densham & Charles 1 198^ 
Codina et al. 1 1984; 
Liu et al. (2006; 
Reig et al. 1 1996' 
Negueruela et al. 1 19961 
Negueruela & Reig I 200l] 
Wolff & Morrison 1 1974] 
Bonnet-Bidaud & Moucnet ||19£ 
Walborn 1 1973! 
Klegueruela 1 1998j 
Liu et al. 12006' 
kee et al. 1 1995' 
^Belczynski & Ziolkowski I 20091 
prampton & Hutchings 1 1978j_ 
Hammerschlag-Hensberge et al. 
l^arkes et al. 1 1980; ' 
Israel et al. ^2000^ 



Negueruela & Okazaki "f200l1 
pkazaki & Negueruela |200i; 
Belczynski fc Ziolkowski J,2009] 
Grampton et al. |(l978j 



Negueruela fc Reig 1 )2004] 
kaur et al. I 2008; 
keig et al. 1 2010! 
karba et al. 1 2006] 
Liu et al. 12006: 

Tomsick et al. 120 0^ 

Torreion & N eg ueruela ||2004] 
Hahoui et ar| i |20D§] 
Goleiro et al. (in pr eparation') 



Butler et al. '('20081^ 

Nespoli et al. 1200 7] 

'Zurita Heras & Walter I J2OO9] 
Negueruela et al. I 2003; 
Johnston et al. 1 1994; 
keig et al. (2005; 
lleig et al. I 200 1; 
Liu et al. (2006: 
Liu et al. 12006; 

keig et al. I 2004^ ^ 

krinia & Massa I 2010; 
Negueruela et al. ( 2008J. 



Belczynski & Ziolkowski (2009; 



[Belczynski 



; Ziolkowski J20091 



TABLE 1 

The HMXB sample. Distance D, Extinction Ay, T, R, spectral type SpT a nd spectral type re f erence SpT Ref . ar e given. 

D AND Ay AR E COMPUTED IN THIS ARTICLE WHEREAS T AND R ARE DERIVED FROM 'VACCA ET AL.I (|1996|) , |PANAGIAI (Il973l) AND [SeARLE ] 
I ET AL.| ( |2Q08[ ) OR TAKEN IN THE LITERATURE WHEN AVAILABLE (SEE TABLE [2]) . WHEN DIFFERENT SPECTRAL TYPES ARE PROPOSED IN THE 
LITERATURE, THE DIFFERENT R AND T VALUES USED ARE GIVEN. SPECTRAL TYPES BETWEEN BRACKETS ARE THE SPEC TRAL TYPES USED 
FOR R AND T DETERMINATION. D AND Ay ERROR DETERMINATION ARE DEVELOPED IN SECTION |2.2| 



SOURCE NAME R (Rp) 



T(K) 



Reference 

"Reig et all (|l996t 
Herrero et aL ( 1995} 
Stee et aL (1995) and [Goraya] ( |2007| > 
Kaper et aL (20061 
Reynolds et al. (199^^ 
Crampton et al. (1978 ) 
Negueruela et al. ( 20li| > 
Reig et al. (2001) 



3A 0114+650 
Cyg X-1 
Gam Cas 
GX 301-2 

H 1538-522 (reL 1) 
H 1538-522 (reL 2) 
PSR B1259-63 
RX J0812.4-3114 



37.0 ± 15.0 
17.0 
10.0 
62.0 

17.2 ± 1 
17 ± 2.0 

10 ± 2 



24000 ± 3000 
32000 

25000-30000 
20400 

28000 ± 2000 
31500 ± 1000 

32000l?ro 
28000 ± 2000 



TABLE 2 

Radius R and Temperature T values available in the literature with references. 
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AT A TV /TTT^ 

NAME 


D (kpc) 


Av (mag) 


Dref (kpc) 


Avref (mag) 


lA 0535+262 


3.8 ± 0.33 


1.9 ± 0.26 


1.8 ± 0.6 


2.7 ± 0.2 


lA 1118-615 


3.2 ± 1.4 


4.6 ± 1.9 


3.6 ± 0.9 


4.2 max 


IE 1145.1-6141 


10.5 ± 0.90 


5.2 ± 0.19 


8.2 ± 1.5 


- 


IH 1249-637 


0.63 ± 2.5 


1.8 ± 2.5 


0.43 ± 0.060 

n Q,r\max=0.36 
^•'^^min=0.26 


- 
- 


3A 0114+650 


6.5 ± 3.0 


4.0 ± 0.60 


7.0 ± 3.6 


- 


3A 0726-260 


5.0 ± 0.82 


2.5 ± 0.25 


6.1 ± 0.30 




3A 2206+543 


3.4 ± 0.35 


1.8 ± 0.60 


2.6 


- 


4U 1700-377 


1.8 ± 0.15 


2.0 ± 0.15 


1.7 




Cep X-4 


3.7 ± 0.52 


5.3 ± 1.4 


3.8 ± 0.60 


- 


Cyg X-1 


1.8 ± 0.56 


3.4 ± 0.18 


-0.11 


- 


EXO 0331+530 


6.9 ± 0.71 


6.0 ± 0.50 


6 < d < 9 


- 


EXO 2030+375 


3.1 ± 0.38 


12 ± 1.4 


7.1 ± 0.20 




gam Cas 


0.17 ± 0.50 


1.2 ± 2.2 


1 SS'"^O'^=0-208 
min=0.168 


- 


GRO J1008-57 


4.1 ± 0.59 


6.7 ± 1.1 


5 


- 


GT 0236+610 


1.8 ± 0.20 


3.8 ± 0.65 


2 


- 


GX 301-2 


3.1 ± 0.64 


6.3 ± 0.14 


3-4 


- 


GX 304-1 


1.3 ± 0.10 


6.0 ± 0.17 


2.4 ± 0.50 


- 


Ginga 0834-430 


7.1 ± 4.2 


11 ± 2.2 


3 < d < 5 


- 


H 0115+634 


5.3 ± 0.44 


6.4 ± 0.28 


7-8 


- 


H 1145-619 


4.3 ± 0.52 


1.7 ± 0.55 


3.1 


- 


H 1417-624 


7.0 ± 0.74 


6.1 ± 1.1 


1.4 < d < 11.1 


6.1 < Av < 8.9 


H 1538-522 


6.2 ± 1.8 


6.4 ± 0.28 


5.5 ± 1.5 
6.4 ± 1.0 
4.5 


- 
- 

6.5 ± 0.3 


IGR J00370+6122 


3.4 ± 1.2 


2.4 ± 0.19 


3.0 


- 


IGR J01583+6713 


4.1 ± 0.63 


4.7 ± 0.32 


4.0 ± 0.4 


4.5 ± 0.2 


IGR J08408-4503 


3.4 ± 0.35 


1.7 ± 0.60 


2.7 


- 


IGR J11215-5952 


7.3 ± 0.68 


2.6 ± 0.60 


8.0 


- 


IGR J11305-6256 


3.6 ± 0.71 


1.0 ± 2.2 


3 


- 


IGR J11435-6109 


9.8 ± 0.86 


5.7 ± 0.28 


4.5 




IGR J16465-4507 


12.7 ± 1.3 


5.0 ± 0.75 


12.5 




IGR J18214-1318 


10.6 ± 5.0 


14 ± 3.1 


10 




IGR J18410-0535 


7.8 ± 0.74 


6.1 ± 0.65 


3.2 




IGR J18450-0435 


6.4 ± 0.76 


6.7 ± 0.49 


3.6 




KS 1947+300 


8.5 ± 2.3 


4.0 ± 0.49 


10 


3.38 


PSR B1259-63 


1.7 ± 0.56 


3.8 ± 0.70 


2.3 




RX J0440.9+4431 


2.9 ± 0.37 


2.9 ± 0.25 


3.3 




RX J0812.4-3114 


8.6 ± 1.8 


2.3 ± 0.20 


8.8 




RX J1744.7-2713 


1.2 ± 0.46 


2.7 ± 0.57 


0.8 ± 0.1 




SAX J1818.6-1703 


2.7 ± 0.28 


1.6 ± 0.80 


2.5 




SAX J2103.5+4545 


8.0 ± 0.78 


4.2 ± 0.25 


6.5 


4.2 ± 0.3 


Vela X-1 


2.2 ± 0.22 


2.2 ± 0.46 


1.9 ± 0.1 




XTE J1855-026 


10.8 ± 1.0 


5.8 ± 0.90 
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XTE J1946+274 


6.2 ± 3.0 


6.9 ± 0.74 


9.5 ± 2.9 




X Per 


1.2 ± 0.16 


0.81 ± 0.22 


0.70 ± 0.30 


1.05 ± 0.02 



D and Ay Ref. 

Giangrande et a E] (|lp8Q|> 



Janot-Pacheco et 'al.Vl98l| 
Uensham & Charles ( 1982| 
Uodina et al. (1984) 
Uhevalier & Ilovaisky^ ( 1998| > 

Reig et al. (1996) 

Negueruela et al. ( 1996| > 
Blay et al. (2006j 
Xnkay et al. (2001 



Bonnet-Bidaud 
Reid et al. (2011)" 
Negueruela et al.l (11 



ouchet 



(1998]) 



1999 



Wilson et al. (20 02r 

Uhevalier & Ilovaisky ( 1998| ) 
Uoe et al. (1994) 
Steele et al. ( 199^ 
gaper et al. (200^) 
Farkes et al. {19^ 
grael et al. (200(F 



Negueruela & Okazaki (2001 ) 
Stevens et al. (1997) 
Grindlay et al. ( 1984) 
fcrampton et al. ( 1978| 
Reynolds et al^ (^1992| 

glark (2004| 

Reig et al. ( 20051 
gaur et al. (20081 
geyder et al. (200?) 
Negueruela et al. ( 200^ 
Masetti et al. (2006) 



Torrejon & Negueru ela| ( |2004 ) 

Hmith (2004j)^ 

Butler et al.1(|2009t 
Nespoli et al. (2008 ) 
Zurita Heras & Walter 
Kegueruela et al 



Reig et al. 


(2005 


Reig et al. 


('2001! 



HTdoli et al 
Feig et al. 



1. (1997| 
, (200f 
(20045; 




Hadakane et al. 
Uorbet et al. (1999) 
Wilson et al. (2003| 
,Lyubimkov et al" ^997 1 



TABLE 3 

Comparison of D and Ay values derived in this study with Dref and Avref available in the literature. 
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D = I I ± 0.90 kpc 
Av = 5.2 ± 0. 1 9 mag 



Fig. 1. — Result of the fitting for the source IE 1145.1-6141 with 
distance D and extinction in V band Ay- 

appears to be less important for supergiant stars than 
for Zero- Age Main Sequence stars. 

Uncertainties on the radius and temperature of the 
companion star could have a severe influence on distance 
and Av error. We generate SEDs of various sources 
using a large range of temperatures which clearly show 
that the distance error is less sensible to the temperature 
error than to the radius one. That can be easily un- 
derstood given the fact that distance is directly derived 
using the R/D ratio. We distinguished different cases. 
1) For sources for which spectral type is accurately 
determined , we select R and T values mostly in [Vacca 
|et al.||1996| for consistency. If no data is given in this ref- 
erence, we used R and T values fflven in Panagia (1973) 



(|^UU8[ ). With the 



Martins fc Plez| ( |2QQ6[ ) or [Searle et al.| 
aim ot deriving accurate distance and Ay errors due 
to radius and temperature uncertainties, we fitted the 
SED of each system using the R and T values given in 
PVMS and we compute the dispersion in distance and 
Av values obtained using these four references. These 
computations lead to a mean difference in distance value 
of 10% as well as for the Av value. Then, for systems 
for which spectral type is well known, we consider a 
systematic error of 10% of the derived distance and Av 
values due to radius and temperature uncertainties. 2) 
For sources for which values of R and T are available 
in the literature (see Table j2|, we use these values in 
our fitting procedure to derive distance and extinction 
in V band. When maximum and minimum radius and 
temperature values are available, we derive distances 
and dispersion on distance using these values. We also 
derive the distance considering R and T values given 
in PVMS. Then, we consider as the final error (due to 
radius and temperature uncertainties) the dispersion on 
all the distance values obtained with R and T values 
given in the four references and with R and T values 
given in the paper dedicated to the source (referenced in 
Table 2 . The same approach was used to derive an Av 
error due to radius and temperature uncertainties. 3) 
When the spectral type is not well constrained and if no 
R and T value is available in the literature, we derived 
distances (and then a mean distance) using the different 
possible spectral types of the companion star and then 



the different R and T values given in PVMS. We finally 
compute the dispersion on all the distances derived and 
considered it as the error due to radius and temperature 
uncertainties. We derive an Av error due to radius and 
temperature uncertainties for these sources, using the 
same method. 4) For luminosity classes II and IV, since 
no data are given in PVMS we used intermediate R and 
T values between class I and III for class II and between 
class III and V for luminosity class IV. Errors were 
computed using the dispersion on distance obtained 
with these R and T values and distances obtained with 
R and T values of lower and higher luminosity classes. 

Degeneracy between several parameter values (based 
on the fitting procedure) needs to be taken into ac- 
count. Indeed, solely relying on a single best fit does 
not capture the full phenomenology associated with 
SED fitting because D and Av are degenerate in this 
approach. In order to produce the best set of fits and 
to determine the dispersion on distance and extinction, 
we carried out 500 Monte Carlo simulations for each 
observed source, by varying the photometry within the 
uncertainties. Hence, we generate a random number 
from a normal distribution (assuming the photometric 
errors to be Gaussian), defined by the error bars for 
each photometric point, so that we build 500 new SEDs 
derived from the original one. These 500 new SEDs 
are subject to the same statistic computation as 
the one described above. Then, we have an entire set 
of best fits of parameters (D, Av) and we are able to 
plot the distribution in the parameter space, showing 
the distribution of properties derived from these Monte- 
Carlo simulations, and especially showing the dispersion 
on distance and extinction for each source (see Figure 
[2|. This dispersion value is taken as the error of the 
ntting procedure, and the median value of dispersion 
on distance determination for all considered HMXBs is 
then 0.75 kpc. 



5.5 



5.3 



5.2 



5.1 



5.0 




5 



10.0 10.2 10.4 10.6 10.8 11.0 11.2 
D (kpc) 



Fig. 2. — Result of Monte-Carlo simulations for the source IE 
1145.1-6141, to estimate uncertainties in D and Ay due to the 
fitting procedure. Colors represent the number of solutions as a 
function of the two parameters (D and Ay) values. 
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There are other sources of uncertainties, particularly 
the infrared excess of Be stars due to their circumstellar 
envelope generating fre e-free radiation. According to 

this excess should not exceed 
J band, 0.15 magnitude in 
However, this 



Dougherty et al. (|1994| 



a mean 0.1 magnitude m 
H band and 0.25 magnitude in K band, 
value corresponds to absolute magnitude and therefore 
the excess can be smaller in apparent magnitude for 
sources located far away and higher for close ones. To 
take this effect into account, an estimate of the distance 
and extinction is needed. Since these two values are 
derived from the fitting procedure, we are only able to 
consider the distance and absorption values obtained 
without taking this IR excess into account. This 
approach is finally equivalent to adding a conservative 
error of 0.1, 0.15 and 0.25 magnitude to the apparent 
magnitude in the J, H and Ks bands respectively; this 
method presenting the advantage of treating all sources 
in an uniform way. Based on these results, presented in 
Figure |3j we assume that this uncertainty will not affect 
the distances derived in this paper. 

Furthermore, to test the potential additional error 
due to a different extinction l a w, w e use Ry values 



given by Geminale & Popowski (2004). We affected to 
each HMXB an Ry value derived from the closest star 
present in their catalogue and we fitted the SED using 
the Cardelli law with the new Ry value. We find a 
median difference in distance determination of 0.07 kpc. 
Since deriving an accurate Ry value for each system 
appears to be observationaly biased (it depends on the 
closest star available in their catalogue), we decided not 
to take into account this additional uncertainty. 

Finally, Table [l] presents the 46 sources and their fun- 
damental parameters (computed in this work or taken in 
the literature). Figure [To] represents most of the studied 
HMXBs with the uncertainties on their location com- 
puted taking into account all the errors described above: 
the error on distance due to uncertainties on the com- 
panion star radius and temperature and the error com- 
ing from the fitting procedure. We derived with the same 
way, the final error on Ay. 

3. RESULTS: HMXB DISTRIBUTION AND CORRELATION 
WITH STAR FORMING COMPLEXES 

We present in Figure the distribution of HMXBs in 
the Galaxy, obtained wrth our novel approach based on 
distance determination. The spiral arms model given 



by Russeil 



arises: is t 



(2003) is also presented. The question then 
Jnere a correlation between this distribution of 
HMXBs and the distribution of Sta r Forming Complexes 
(SFCs) in the Milky Way (given b y |Russeil|2003 ), as it is 
expected from the short HMXB lifetime Y 'I'he first ap- 
proach we adopt is to carry out a Kolmogorov-Smirnov 
test (KS-test) on each axis in order to quantify whether 
the two samples are drawn from the same probability 
distribution. We obtain a value of 0.15 for the X axis, 
a value of 0.25 for the Y axis and a value of 0.31 for 
the galactic longitude. These values are not negligible, 
suggesting that a correlation between the two samples 
does exist, though part of the information is lost because 
of the projection on the two axis. To overcome this 
caveat, we propose another method described hereafter. 



m 



-15 -10 -5 5 10 15 

X (kpc) 

Fig. 3. — Green circles represent the initial positions of Be stars 
whereas blue stars represent the source positions taking into ac- 
count the IR excess. Red star at (8.5 ; 0) represents the Sun 
location. Dotted pink circles represent radius of respectively 3, 6, 
9, 12 and 15 kpc from the Galactic center. 



Carina 









Criix: , 






X (kpc) 



Fig. 4. — Distribution of HMXBs (blue stars) and SFCs (green 
circles). The circle radius of SFCs represents the excitation pa- 
rameter value. The spiral arm model from Russeil (2003) is also 
plotted and the red star at (8.5 ; 0) represents the Sun position. 

We assume that each HMXB (blue stars on Figure [5| 
is clustered with several SFCs (green circles). Hence, we 
can define two characteristic scales: a typical cluster size 
and a typical distance between clusters. Around each 
HMXB, we define several circles with different radii (red 
circles) and we finally count the number of HMXBs for 
which at least one SFC is within the specified radius 
(called 'number of correlations" hereafter). 

The number of correlations versus the circle radius is 
plotted as the blue curve on Figure |6] The green curve 
is that expected from chance correlations assuming the 
HMXBs are evenly distributed across the sky following 
an uniform distribution on x and y axis between and 15 
kpc. The dashed line represents the difference between 
the two previous curves. This difference, being non equal 
to zero, allows us to state that a strong correlation exists 
between HMXB and SFC positions in the Milky Way. 
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Fig. 5. — Description of the method used to evaluate the corre- 
lation between HMXBs and SFCs. 



HMXBs 



Random 



Difference 




10"^ 10"^ 0.3 10° 1 .7 10^ 10^ 

radius (kpc) 

Fig. 6. — Result of the correlation determination in 2D. 



Moreover, we compute the two characteristic scales de- 
scribed above: the typical cluster size of 0.3 kpc and the 
typical distance between clusters of 1.7 kpc with uncer- 
tainties of 0.05 kpc and 0.3 kpc respectively, using the 
uncertainty on HMXB distance of 17% as derived in Sec- 
tion O 

If we take into account the uncertainties in HMXB po- 
sitions (described in Section 2.2) and in SFC positions 
(given in Russeil 2003, median error of 0.25 kpc), the 
correlation still exists with the s ame cluster size and the 
same distance between clusters. iBodaghee et al. (2012) 



mention that HMXBs and OB complexes are clustered 
for a cluster size r < 1 kpc. This upper limit, obtained 
with a different method, is consistent with our results. 
Finally, we test o ur correlation cod e using a sample of 
globular clusters ( |Bica et al. [ 2006 ), principally located 
in the Galactic bulge, l^lgure [7| shows the result of cor- 
relation test. The number of correlation as a function 
of the circle radius (blue curve) follows the trend of the 
green curve, showing the evolution expected from chance 



10" 
radius (kpc) 



Fig. 7. 
clusters. 



Result of the correlation determination with globular 



4. IMPLICATION OF THE CORRELATION ON HMXB 
FORMATION AND EVOLUTION 

The distribution of HMXBs reflects the stellar forma- 
tion that took place some tens of Myr ago, since they 
are not an instantaneous star formati on rate (SFR) 
indicator as explained in Shtykovskiy & Gilfanov|(|2 007[ ). 
Then, an offset between spiral arms (an indicator of 
the actual star formation) and HMXB distribution is 
expected. Indeed, since the spiral arm rotation velocity 
is different than the angular velocity of the stellar 
disk, we expect HMXB positions to be offset from the 
cu rrently visible SFCs. T his t ime lag was menti oned 



Lutovinov et al. (2005) and Dean et al. (2005) but 
a deeper investigation oi this issue was not possible 
due to the limited samp l e of H MXBs at that time. 
Shtykovskiy fc Gilfanov] (2007) evaluated the offset 
lor the galaxy M51. Here, we attempt to implement 
this formalism in the case of the Milky Way with th e 
approach we have mentioned in Coleiro & Chaty (2011 ) . 



One can suppose the spiral arm density wave t o rotate 



Dias & 



with the speed Qp = 24 kms~^kpc~^ (see e.g. 
[Lupine 2005), in the same direction as the stellar disk, 
which velocity curve is assumed to be flat in the ga lac- 
tocentric distance range of interest, according to [Brand 
fc Blitz] ( |T993l ): 



l//(220km.s"^) = ai(r/8.5kpc)^2 ^ 



(3) 



with ai=1.00767, a2=0.0394, a3=0.00712 and r is the 
radius from the Galactic center. 

Then, to a first approximation, it is possible to locate 
the expected HMXB locations relative to the current po- 
sition of the spiral arms in time r, i.e., the angular offset 
A9(r) by the equation: 



A6(r) = {n{r)-np)r, 



(4) 
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where Q(r) is the galactic rotation curve derived from 
Equation [3j To estimate the displacement of HMXBs 
relative to the current position of the spiral density 
wave, we plot the angular offset, A6(r), as a function of 
the radius for sources with different ages (10, 20, 40, 60, 
80 and 100 Myr.) (Figure |8|. In a second time, we plot 
the expected positions of sources formed 10, 20, 40 Myr 
ago in the Galaxy map (Figure [9| (For better visibility, 
here, we do not plot the expected positions of sources 
formed 60, 80 and 100 Myr ago). 



10 Myr 
20 Myr 
40 Myr 
60 Myr 
80 Myr 
100 Myr 




Fig. 8. — Evolution of the offset angle between spiral arms and 
expected positions of sources as a function of the radius in the 
Galaxy according to the age of sources. The vertical dashed line 
indicates the corotation radius. 



Even if Figures [9] and 10 do not visually suggest any 
substantial offset between the current spiral arm position 
and the expected position of HMXBs (which depends on 
the age of the sources) we would now like to quantita- 
tively assess the presence of an offset. 

4.1. Existence of an offset between HMXBs and 
Galactic spiral arms 

To perform this study we need to consider a number 
of different issues. As underlined above, the expected 
offset depends on the age of the X-ray sources, hence, to 
highlight this offset, we must split the sample of HMXBs 
depending on the age of the sources. Two different sam- 
ples are then created: one containing 4 supergiant stars 
(luminosity classes I or II according to Charles & Coe 
12006) and a second one contai ning 9 Be stars (lu minos- 
ity class HI or V according to [Charles fc Coe|[2Q06p . We 
explain the way these samples were created in Section 
K2\ . 

We calculate the distance from each HMXB to the clos- 
est actual spiral arm (given in Russeil 2003) using: 



distance = \/ (Xg 



)2 + (Y 

S( 



source 7 -^source 
'-arm? ^rm 



^ arm / i 

(5) 

) are the coordinates of the 
) are the coordinates of the closest 



where (Xg, 
source and (Xa 
point on the arm. 

We follow exactly the same procedure to calculate the 
distance from each source to the closest expected position 




Fig. 9. — Distribution of HMXBs and spiral structure. The blue, 
green and red curves correspond respectively to the expected lo- 
cations of objects with ages of 10, 20 and 40 Myr computed with 
Equation pi For better clarity, we do not represent expected loca- 
tions of objects with ages of 60, 80 and 100 Myr. The grey curve 
represents the current position of spiral arms and the dashed circle 
indicates the corotation radius. Dotted circles represent radii of 3, 
6, 9, 1 2 an d 15 kpc. Dashed frame represents the region shown in 
Figure [To] 

of sources formed 20, 40, 60, 80 and 100 Myr ago and 
we determine the mean value of the offsets (taking into 
account all the sources of the two samples) . The sources 
are expected to be closest to one of the expected positions 
computed above than to the current spiral arms observed 



by Russeil (2003) . The method is described in Figure 



(for mstance, m Figures 10 and 11, a source of 20 ]V 



n 

yr 

should be located closer to the green arm representing 
the expected position of a 20 Myr-old HMXB than to 
the current spiral wave position). Results are given in 
Figure 12] 

We determine the error bars using the propagation of 
uncertainty formula. The 1-a error associated to the dis- 
tance between each source i and the closest point on the 
arm is then given by the following equation: 



2 

^source 



^ V 



{Xs. 



{Ysc 



arm ) 
2 



■ (7 



(6) 



-^arm)^ H~ (^c 



-'- arm / 



where {ax^, 



are the errors on source po- 
sition detailed in Section |2. 2 and {o'x^^^ , ^y^^^) are the 



O'Y 



errors on arm position taken to be equal to zero. Finally, 
the uncertainty associated to the mean value is calculated 
as follows: 



N 



(7) 



where N is the number of sources and ai is the error 
associated with the i*^ source. 



ure 



Taking into account all the sources of the sample (Fig- 
lower panel), we cannot observe any significant 



12 



variation of the offset with time before 60 Myr. This re- 
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X (kpc) 

Fig. 10. — Positions of HMXBs with error bars (zoom in the Solar region defined by dashed frame in Figure 19^ 




> T (Myr) 



Fig. 11. — Method used to derive the age of a sample of HMXBs. As in Figure a spiral arm is represented with the expected positions 
of sources of 10, 20 and 40 Myr old. An HMXB is represented by the star and distances to current spiral arm and to each expected positions 
are called di, d2, ds and d4 respectively. The evolution of this distance as a fonction of the time is represented in the diagram on the right. 
In this case, the source would be ~ 20 Myr old. 
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suit was expected because the minimum peak of distance 
between each HMXB and the closest expected position 
should be clearly different whether we consider Be or su- 
pergiant stars. Also, by taking all spectral types into ac- 
count, we tend to loose a part of the information except 
the mean age upper limit of 60 Myr highlighted by the 
plot. For both super giant and Be samples, we observe a 
more significant increase of the offset after 60 Myr (Fig- 



ure [12, first and mid panels). This increase could be a 
signature of the expected offset between the HMXB po- 
sitions and the current spiral arms. However, we must be 
cautious about this result especially because of the small 
number of HMXBs in the two samples. 

Even if the locations of the sources are accurately 
determined, several reasons may affect the HMXB 
de nsity and prevent the o ffset detection as underlined 
by Lutovinov et al. ([2005): complex motion of density 
wave and stars froni their birthdate to the X-ray phase, 
presence of previously undetected parts of the Galactic 
spiral arms, observational selection effect, etc. Finally, 
a larger sample of supergiant type HMXBs is needed to 
confirm the offset detection more confidently. 

Moreover, the time interval during which HMXBs ap- 
pear should transla te the mass rang e of both stars of 
binary systems (see Dean et al. 2005). The results pre- 
sented on Figure p!2] only enable us to state that this time 
interval is lower than 60 Myr on average for all stars. 

4.2. Deriving the age of HMXBs 

Following the method described above, one can 
compute a distance between each source and the dif- 
ferent theoretical arms that correspond to the current 
predicted position of a source sample born 20, 40, 60, 
80 or 100 Myr ago. In the previous section we studied 
a sample of Be and supergiant HMXBs, and this could 
also be applied to each source separately, to constrain 
the age and the potential migration of the system due 
to a supernova kick. 

We choose not to take into account here the arm 
width that could be translated as an uncertainty on the 
position of the expected position of sources. Indeed, by 
taking this position dispersion into consideration, 1-a 
errors considerably increase and prevent any conclusion. 
Then, we assume all sources to be formed at the central 
position of the arm. 



it 
O 



Be 




T (Myr) 

Fig. 12. — Evolution of the mean distance (in kpc) between source 
and closest expected position versus time r (in Myr) for Be sources 
(upper panel), supergiant sources (mid panel) and all the sources 
(lower panel). Error bars are 1-a uncertainties. 

cal positions observed for some sources. Since off- 
sets between expected positions of sources of differ- 
ent ages are very small, it can be possible to ascribe 
an incorrect age to a source located close to an ex- 
pected position that does not correspond to its real 
age; 

3. sources that underwent more significant kicks could 
be incorrectly associated with a theoretical position 
of sources of different age. 

Then, for this study, we only take into account sources 
presenting an expected offset evolution during time, to 
make the conclusions easier. 



We present on Figure 13 results for BeHMXBs and 
sgHMXBs for which the evolution of the offset as a 
function of time is consistent with the theory. It is 
possible to determine the rough age of these sources 
and a lower and upper limit of kick migration distances. 
The results are summarized in Table [4] and will be 
commented in the following paragraphs. 



We expect the distance from the source to the expected 
position to decrease until the expected position corre- 
sponds to the age of the source, and then to increase 
afterwards (see Figure 11). For some sources, the off- 
set change is more complicated and can be explained as 
follows: 

1. belonging to one of the four arms of the Milky Way 
is not well established when the distance between 
a source and the theoretical expected positions is 
seen as always increasing with time; 

2. for sources located close to the corotation ra- 
dius, current spiral arms and expected positions of 
HMXBs of different ages are almost superimposed 
explaining the quasi- const ant distance to theoreti- 



If we consider only the sources studied separately here 
(sources showing the expected evolution of offset as a 
function of age), we can derive a mean age of ~ 45 Myr 
for SgHMXBs and 51 Myr for BeHMXBs. This result 
seems consistent with the distinct evolution timescales 
of these two kinds of systems but again, we must be 
cautious about this result because of the small sample of 
sources used in this calculation (4 supergiant and 9 Be 
systems). 

4.3. Constraints on migration from supernova kicks 

Even if the theory is still poorly understood, it is 
admitted that a small asymmetry in the way a super- 
nova explodes could make a binary system move due 
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SOURCE NAME AGE (Myr) MIGRATION DISTANCE (kpc) UNCERTAINTY 







Be 






lA 0535+262 


80 




0.10 


0.30 


lA 1118-615 


80 




0.088 


0.56 


EXO 0331+530 


60 




0.25 


0.080 


GRO J1008-57 


40 




0.074 


0.15 


GX 304-1 


40 




0.048 


0.59 


H 1417-624 


20 




0.20 


0.39 


PSR B1259-63 


60 




0.037 


0.51 


RX J0440. 9+4431 


20 




0.011 


0.17 


RX J1744.7-2713 


60 




0.10 


1.0 






Supergiants 






4U 1700-377 


80 




0.15 


0.28 


ICR J16465-4507 


20 




0.087 


0.052 


ICR J18410-0535 


60 




0.013 


0.11 


H 1538-522 


20 




0.14 


0.52 



TABLE 4 

Age and migration distance derived respectively for 
BeHMXBs (top) and sgHMXBs (bottom). Uncertainties 
are derived by equation [6] i.e. correspond in figure [it] to 
the error on the minimum distance d3 . 



to a substantial increase of its velocity in a particular 
direction. We computed the cluster size (in the same 
way as in Section [3j only for Be stars that are expected 
to have experienced a supernova kick event during their 
evolution and we found a value of 0.3 kpc, consistent 
with the migration distances underlined in Bodaghee 
|et al. (2012). This value also enables us to constrain the 
time lapse between supernova event and HMXB stage, 
assuming a regular value of kick velocity of 100 kms~^ 
(for hi gher kick v elocity, fewer system should remain 
bound. Hills 1983) and a maximum migration distance 
due to the kick event of 0.3 kpc, we reach an upper limit 
for the time lapse between the supernova explosion and 
the HMXB step of around 3 Myr. 

To improve this study we focus on the sources selected 
before (4 sgHMXBs and 9 BeHMXBs, Table [4|. The 
distance between the object and the closest expected po- 
sition gives a kick value for each source and a mean value 
depending on the spectral type of the two samples. Re- 
sults are presented in Table [4j We derive a mean migra- 
tion distance of 0.11 kpc for BeHMXBs and of 0.10 kpc 
for SgHMXBs. Again, we should be careful about these 
mean values because of large error bars and small sam- 
ples. Moreover, it is important to underline that these 
derived values only represent a lower limit to the kick 
migration distance since we cannot take into account the 
migration distance on galactic latitude given by the kick. 
Then, we only get the projected migration distance on 
the Galactic plane. 

5. CONCLUSION 

Examining the distribution of HMXBs is of major in- 
terest in order to study in depth the formation of these 
high energy sources. However, HMXB locations are usu- 
ally poorly constrained and largely dependent on the de- 
termination method. Here, for the first time, we de- 
termine the location of a sample of HMXBs using an 



uniform and accurate approach: SED fitting of their dis- 
tance and absorption. This method, based on a least- 
square minimization, enables us to reveal a consistent 
picture of the HMXB distribution, following the spiral 
arm structure of the Galaxy. The uncertainties lead to 
a small error on source location and allow us to tackle 
the study of the correlation with SEC distribution. This 
study shows that HMXBs are clustered with SFCs and 
enables to quantitatively define the cluster size (0.3 ± 
0.05 kpc) and the distance between clusters (1.7 ± 0.3 
kpc). We go further by quantitatively assessing the off- 
set between current spiral density wave position and ex- 
pected HMXB positions due to the fact that the matter 
rotation velocity is different than the spiral arm rota- 
tion speed. Exploring the environment in which such 
binary systems were formed is of major interest to study 
the properties of these binary systems such as stellar 
mass, dust cocoon density, etc. Here, we quantitatively 
show the correlation between HMXB distribution and 
Star Forming Complexes distribution. Even if we high- 
lighted the expected offset between current spiral arms 
and source positions for some sources, it remains dif- 
ficult to assess it for the entire sample. Undoubtedly, 
this assessment will be improved using a larger sample 
of sources, accurate Galactic spiral arms model and a dy- 
namical model of matter and density wave. Our method 
of investigation does not give exhaustive results for the 
entire sample of sources because several sources are lo- 
cated close to the corotation radius and, for some sources, 
an association with one of the four arms of the Milky Way 
is not well established. However, for 4 sgHMXBs and 9 
BeHMXBs, we are able to derive an age (mean age of 51 
Myr for BeHMXBs and 45 Myr for sgHMXBs) and a mi- 
gration distance (mean value of 0.11 kpc for BeHMXBs 
and 0.10 kpc for sgHMXBs), giving constraints on the 
supernova explosion kick. This study represents an im- 
portant progress in the investigation of formation and 
evolution of these binary systems. 
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